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Abstract 

Consistency of Einstein-de Sitter model with the recently observed 
SNe la by the Hubble Space Telescope is examined. The model shows 
a reasonable fit to the observation, if one takes into account the ex- 
tinction of SNe light by the intergalactic metallic dust ejected from 
the SNe explosions. Though the fit to the new data is worsened con- 
siderably compared with the earlier data, it can still be regarded ac- 
ceptable. We should wait for more accurate observations at higher 
redshifts (as expected from the coming space missions such as SNAP 
and JWST) in order to rule out a model, which seems to explain all the 
other existing observations well (some even better than the favoured 
ACDM model), is consistent with beautiful theoretical ideas like in- 
flation and cold dark matter, and is not as speculative as the models 
of dark energy. 

Subject heading: cosmology: theory, whisker dust, SNe la. 
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1. INTRODUCTION 

The high redshift supernovae (SNe) la explosions look fainter than they are 
expected in the Einstein-de Sitter (EdS) model (flat cold dark matter model 
with fitotai — 1)) which used to be the favoured model before these obser- 
vations were made a few years ago. The most popular explanation of this 
observed faintness is given by invoking 'dark energy', a homogeneously dis- 
tributed energy component with negative pressure. This happens because 
the metric distance of an object, out to any redshift, can be increased by 
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incorporating a 'fluid' with negative pressure in Einstein's equations. The 
simplest and the most favoured candidate of dark energy is Einstein's cos- 
mological constant A. 

Alternative mechanisms have also been suggested in order to explain the 
observed faintness of the distant SNe la, without invoking cosmic acceler- 
ation. One of such attempts considers photon-axion oscillations in the in- 
tergalactic magnetic field as a way of rendering supernovae dimmer without 
invoking dark energy (Csaki, et al 2002; 2004). 

Another alternative way to explain the faintness of the high redshift SNe 
la is to consider the absorption of light by metallic dust ejected from the SNe 
explosions— an issue which is generally ignored while discussing m-z relation 
for SNe la. It was shown earlier (Vishwakarma, 2003) that if this effect is 
taken into account, the EdS model explains the SNe la data from Perlmutter, 
et al (1999) (together with SN 1997ff at z « 1.7) successfully. Since then, 
many SNe la at higher redshifts have been discovered with the Hubble Space 
Telescope (HST) (Riess, et al 2004). Extending our earlier work further, we 
examine in this paper how well or badly the EdS model fits the new data 
which include 7 highest redshift SNe la known so far, all at z > 1.25. 

2. EXTINCTION BY METALLIC DUST 

Chitre and Narlikar (1976) were the first to discuss the role of intergalactic 
dust in the m-z relation, which was however largely ignored that time. It is, 
however, generally accepted now that the metallic vapours are ejected from 
the SNe explosions which are subsequently pushed out of the galaxy through 
pressure of shock waves (Hoyle & Wickramasinghe, 1988; Narlikar et al, 
1997). Experiments have shown that metallic vapours on cooling, condense 
into elongated whiskers of ks 0.5 — 1 mm length and ^10~ 6 cm cross-sectional 
radius (Donn & Sears, 1963; Nabarro & Jackson, 1958). It has been shown, 
by considering standard big bang and alternative cosmologies, that this kind 
of dust extinguishes radiation travelling over long distances (Aguire, 1999; 
Banerjee, et al 2000; Narlikar, et al 2002; Vishwakarma, 2002; 2003). The 
dust presumably re-radiates the absorbed radiation at another wavelength. 
For example, it degrades a radiation of high-frequency into the far infrared. 
The imergent radiation is then redshifted into the millimeter waveband due 
to the expansion of the universe. However, this is not the subject of the 
present paper. The density of the dust can be estimated along the lines of 
Hoyle et al (2000): If the metallic whisker production is taken as 0.1 M 
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per SN and if the SN production rate is taken as 1 per 30 years per galaxy, 
the total production per galaxy (of spatial density ~ 1 per 10 75 cm 3 ) in 
10 10 years is p=j 2/3 x 10 41 g. The expected whisker density, hence, becomes 
2/3 x 10 41 x 10~ 75 10~ 34 g cm -3 . We shall see later that this value is in 
striking agreement with the best-fitting value coming from the SNe la data. 

In the following we calculate the mean extinction of the SNe light from 
the whisker dust. Suppose that the light we observe today was emitted at 
the epoch of redshift z from a SN of luminosity L. In traversing a distance 
d£ through the intergalactic medium, the light looses its brightness by an 
amount dL due to the absorption by the whisker grains of mean density p g . 
This is given by 

dL = -k p g L d£, (1) 

where k is the mass absorption coefficient, which is effectively constant over 
a wide range of wavelengths and is of the order 10 5 cm 2 g^ 1 (Wickramasinghe 
& Wallis, 1996). Hence the total loss in the luminosity AL of SN light, in 
traversing a proper distance £(z), can be obtained by integrating equation 

(1): 
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For the Robertson- Walker (RW) metric, this reduces to 

dz' 



AL(z) = exp 
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H(z' 
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where p go is the whisker grain density at the present epoch and p g S 3 = 
Pgo Sq = constant, S(t) being the scale factor of the RW metric. Here and 
henceforth we have considered the speed of light c = 1. The corresponding 
increase in the magnitude Am (= — 2.51ogAL) is thus given by 

fZ Ay 1 

Ara(z) = 1.0857 x k p g0 / (1 + z')° 
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The net magnitude is then given by 
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net 



(z) = m{z) + Am(z), 



(5) 



where the first term on the right corresponds to the usual apparent magnitude 
of the object resulting from the cosmological evolution: 



m(z) =51og[# d L (z)]+M, 



(6) 
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with the luminosity distance di, given by 



for the k — case of the RW metric. The constant M. appearing in equation 
(6) is given by M. = M — 5 logi^o + constant, where M is the absolute 
magnitude of the SNe. For the EdS model (fi m o — 1> ^ao = 0), equation (5) 
reduces to 

m nc \z) = 51og[2{(l+^)-(l+^) 1/2 }]+0.7238x « Pg0 H^l+zf' 2 -1]+M. 

(8) 

3. DATA FITTING 

We consider the data recently published by Riess et al (2004) which, in 
addition to having previously observed SNe, also include 16 newly discovered 
SNe la, 6 of them being among the 7 highest redshift SNe la known, all at 
redshift > 1.25. We particularly focus on their 'gold sample' of 157 SNe la 
which is believed to have a 'high confidence' quality of the spectroscopic and 
photometric record for individual supernovae. We note that the data points 
of this sample are given in terms of distance modulus \x a = m nct — M = 
51ogG?L + constant. However, the zero-point absolute magnitude were set 
arbitrarily for this sample. Therefore, while fitting the data, we can compare 
the observed /i with our predicted m net given by equation (8) and compute 
X 2 from 

m nct (zj; KpgpHp 1 , M)-fJ. ,i 2 ^ 

On - 

The constant M thus plays the role of the normalization constant. The 
quantity o^ o . is the uncertainty in the distance modulus /i 0i j of the i-th SN. 

There are only two free parameters in this model which are to be esti- 
mated from the data: Kp^H^ 1 and M.. The parameter np^H^ 1 , which is 
dimensionless, is of the order of unity if one considers k of the order 10 5 cm 2 
g~\ P g o of the order 10" 34 g cm -3 and H ~ 70 km s _1 Mpc -1 . However, we 
have kept it as a free parameter to be estimated from the data. 

By varying the free parameters of the model, we find that the minimum 
X 2 is obtained for M. = 43.41 and kp^Hq 1 = 4.77 (where k, p g o and H have 
been measured in units of 10 5 cm 2 g _1 , 10~ 34 g cm -3 and 100 km s -1 Mpc -1 
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respectively). These give a value x 2 — 200.99 at 155 degrees of freedom (dof), 
that is, x 2 /dof = 1.29, which represents a reasonable fit. In the absence of 
the extinction from the whisker dust, i.e. for p g0 = 0, one obtains a too high 
X 2 /dof = 324.70/156 = 2.08, with M = 43.58. 

Though there is not a clearly defined value of x 2 /dof for an acceptable fit, 
a 'rule of thumb ' for a moderately good fit is that x 2 should be roughly equal 
to the number of dof. A more quantitative measure for the goodness- of- fit 
is given by the x 2 -P r °bability which is very often met with in the literature 
and its compliment is usually known as the significance level (should not be 
confused with the confidence regions). If the fitted model provides a typical 
value of x 2 as a; at n dof, this probability is given by 



Roughly speaking, it measures the probability that the model does describe 
the data and any discrepancies are mere fluctuations which could have arisen 
by chance. To be more precise, Q(x,n) gives the probability that a model 
which does fit the data at n dof, would give a value of x 2 as large or larger than 
x. If Q is very small, the apparent discrepancies are unlikely to be chance 
fluctuations and the model is ruled out. It may however be noted that the x 2 - 
probability strictly holds only when the models are linear in their parameters 
and the measurement errors are normally distributed. It is though common, 
and usually not too wrong, to assume that the x 2 - distribution holds even 
for models which are not strictly linear in their parameters, and for this 
reason, the models with a probability as low as Q > 0.001 are usually deemed 
acceptable (Press et al, 1986). Models with vastly smaller values of Q, say, 
10~ 18 are rejected. 

The probability Q for the best-fitting EdS model is obtained as 0.008, 
which is though very small, but acceptable. In the absence of the whisker 
dust, this probability reduces to 1.5 x 10~ 12 , which is too small to be accepted. 
In order to compare these results with those in the ACDM cosmology, we 
note that for a constant A, the best-fitting models are obtained as 



Global best-fitting model: fi m0 = 0.46, fi A o = 0.98, M = 43.32, with x 2 /dof 
= 175.04/154 = 1.14, Q = 0.118; 

Best-fitting flat model: fi m0 = 1 - fi A o = 0.31, M = 43.34, with x 2 /dof 
= 177.07/155 = 1.14, Q = 0.109; 




(10) 
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Concordance model (Q m0 = 1 - Ovo = 0.27): x 2 /dof = 178.17/155 = 1.15 
with Q = 0.098. 

We note that the fit to the EdS model, even with the whisker dust, is con- 
siderably worse than those in the ACDM models. However we also note that 
the new data has a general tendency to worsen the fit to any model. For 
example, the earlier data with 57 points (Vishwakarma, 2003) gave better 
fits to the models mentioned above: 

EdS model (with whisker dust): x 2 /dof = 68.97/55 = 1.25, Q = 0.098; 
Best-fitting ACDM: fi m0 = 0.64, fi A0 = 1.20 with x 2 /dof = 57.78/54 = 1.07, 
Q = 0.337; 

Best-fitting flat ACDM: fi m0 = 1 - fi A0 = 0.32 with x 2 /dof = 59.67/55 = 
1.08, Q = 0.31. 

We note that the new data have worsened the fit to the EdS model to a greater 
extent compared to the ACDM model, however, considering the successes of 
the EdS model, we should be tolerant of the low probability associated with 
it, until it goes down to some value like the one obtained without the whisker 
dust. The reasons are many: 

(i) The EdS model is fully consistent with the recent CMB observations made 
by the WMAP. In fact, there is a degeneracy in the f2 m0 — ^ao plane along a 
line f2 m o+f2Ao ~ 1 and a wide range of f2 m o is consistent with the observations 
(Vishwakarma, 2003). For example, with Q^q = 0.1 and Hq = 55 km s _1 
Mpc -1 , this model yields the positions of the first three acoustic peaks at the 
multipole values ipca^ = 220.4, £ P cak 2 — 521.4, £ P cak 3 = 784.9, which are in 
good agreement with WMAP. Additionally a lower value of H is favourable 
for this model to have sufficient age. For example, H should be < 54 km 
s _1 Mpc -1 for EdS model to have an age of the universe > 12 Gyr, so that 
the age of the oldest objects detected so far, e.g., the globular clusters of 
age t GC = 12.5 ± 1.2 Gyr (Cayrel et al. 2001; Gnedin et al. 2001), can 
be accommodated. There are many observations, like the HST observation 
giving H = 0.72 ± 3(stat)±7(systematic) km s _1 Mpc -1 (Freedman et al. 
2001), which give a rather higher value of H . However, there are several 
observations which also measure smaller values of H . For example, there is 
another HST Key Project which gives H = 64±l km s" 1 Mpc" 1 (Jha et al. 
1999). Sandage and his collaborators find a value even as low as H = 58 ±6 
km s -1 Mpc" 1 from an analysis of SNe la distances (Parodi et al. 2000). 
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Additionally, unlike the A-dominated models, the EdS model has no 
strong integrated Sachs- Wolfe effect, so is in better agreement with the low 
quadrupole seen by WMAP (Blanchard, 2003). 

(ii) The recent data on distant x-ray clusters obtained from XMM and Chan- 
dra projects indicate that the observed abundances of clusters at high red- 
shift, taken at face value, give 0.9 < Q m0 < 1.07 (at 1 a) (Blanchard, 2005). 
This is in striking agreement with a matter-dominated EdS model and is 
hard to reconcile with the concordance ACDM model. 

(hi) The estimated value of the parameter kp^Hq 1 from the SNe fit is indeed 
of order of unity, as expected from the theoretical reasoning described in 
section 2. 

(iv) Unlike the models of dark energy, the EdS model is not very speculative. 

In order to have a visual comparison of the fits of different models to 
the actual data points, we magnify their differences by plotting the relative 
magnitude with respect to a fiducial model (fi m0 = 0, f2 A0 = 0, without 
whiskers), [which has a reasonably good fit (% 2 /dof = 191.701/156 = 1.2) 
with M. = 43.40]. This has been shown in the 'modified' Hubble diagram in 
Figure 1. In Figure 2, we have shown the allowed regions in the parameter 
space KpgoH^ 1 — M. at 95% and 99% confidence levels. 



4. BAYESIAN APPROACH 

The frequentist's goodness-of-fit test of models, described in the previous sec- 
tion, uses the best-fitting parameter values to evaluate relative merits of the 
models under consideration. Thus it judges only the maximum likely perfor- 
mance of the models. The Bayesian theory, on the contrary, does not hinge 
upon the best-fitting parameter values and evaluates the overall performance 
of the models. It has thus appeared as a powerful tool of model comparison. 
The Bayes factor, which is a ratio of average likelihoods (rather than the 
maximum likelihoods used for model comparison in frequentist statistics) is 
given by 

y CM) ~ pMMj)' [ } 

where the likelihood for the model Mj, C(M{) is the probability p(D\Mi) to 
obtain the data D if the model Mj is the true one (for more details on the 
Bayesian theory, see Drell et al, 2000; John and Narlikar, 2002). For a model 
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Figure 1: Hubble diagram of the 'gold sample' of SNe la minus a fiducial model 
(f2 m o = 0, ^Ao = 0) : The relative magnitude(Am net = m net — mgduciai) is plotted 
for some best-fitting models. The solid curve corresponds to the EdS model with 
the whisker dust, the dotted curve corresponds to the flat ACDM model, the 
dashed curve corresponds to the spherical ACDM model, and the dashed-dotted 
curve corresponds to the EdS model without any whisker dust. Some SNe with 
redshift between and 1 (mostly with z ~ 0.5) are missed by all the models and 
seem to be outliers. 

Mi with free parameter, say, a and f3 (generalization for the models with 
more parameters is straight forward), this probability is given by 



where p(a\Mi) and p(f3\Mj) are the prior probabilities for the parameters 
a and (3 respectively, assuming that the model Mj is true. £j(a,/3) is the 
likelihood for a and (3 in the model Mj and is usually given by the x 2 -statistic: 



£(Mi)=p(D\Mi) = J da J dp p{a\Mi) p(/3|M,) £,(«,/?), (12) 




(13) 
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Figure 2: The allowed regions by the 'gold sample' of SNe la data (with 157 
points) are shown in the parameter space KpgoH^ 1 — M. of the Einstein de Sitter 
model. The parameters k, p g o and Hq are measured in units of 10 5 cm 2 g _1 , 10~ 
g cm -3 and 100 km s~ x Mpc -1 respectively. The inner ellipse denotes the 95% 
confidence region and the outer one, the 99% confidence region. 



where xl is the x 2 -statistic for the model Mj, like the one for the EdS model 
given by equation (9). For flat prior probabilities for the parameters a and 
/3, i.e., assuming that we have no prior information regarding a and f3 except 
that they lie in some range [a, a + Act] and [/3, (3 + A/3], we have p(a|Mj) = 
1/Aa and p(/3\Mi) = l/A/3. Hence the expression for the likelihood of the 
model Mj reduces to 



C{Mi) 



1 



Act A(3 



pa+Aa /-/3+A/3 




/ / exp 


2 


la J (3 



d/3 dct. 



(14) 



The Bayes factor By, given by (11), which measures the relative merits of 
model Mj over model Mj, is interpreted as follows. If 1 < B^ < 3, there is 
an evidence against Mj when compared with Mj, but it is not worth more 
than a bare mention. If 3 < < 20, the evidence against Mj is definite but 
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not strong. For 20 < Bij < 150, this evidence is strong and for B^ > 150, it 
is very strong. 

In order to compare the A^O cosmology with EdS one, we presume that 
the universe is flat spatially. This gives a fair comparison as well, since both 
the models are similar and have the same number (two) of free parameters. 
As we do not have any prior information on the common parameter Ai, we 
take help from it's best-fitting values, which lie between 43 and 44. Hence we 
assign a flat prior on M that it lies in the range [41, 45]. (This probability 
does not have any consequences though, as it cancels out in the Bayes factor.) 
Similarly, for the other parameter kp^Hq 1 in the EdS model (the best-fitting 
value ~ 4.8), we assume that Kp^H^ 1 e [0, 10]. For the parameter f2 m0 in 
the flat ACDM model (the best-fitting value ~ 0.3), we assume that fl m0 e [0, 
1]. When calculated for the 'gold sample', this gives a Bayes factor favouring 
a model with A ^ over one with A = (EdS) as B = 1.98, which though 
indicates an evidence against the EdS model, but not more than a bare 
mention. 

However, if no prior assumptions are made about the spatial geometry of 
the universe, there is a strong evidence against the EdS model: In this case, 
we assign the prior probabilities f2 m o G [0, 3], ^ao £ [—3, 3], as have also been 
considered by Drell et al (2000) and John & Narlikar (2002). This is because 
the best-fitting values are towards higher side. This gives a Bayes factor for 
ACDM model compared to EdS model as B = 113.23, for a conservative 
prior probability Kp^H^ 1 e [4,5]. This shows a strong evidence against the 
EdS model, though it is not very strong. For a more liberal prior, however, 
this evidence moves towards very strong side. 

One should also note that a proper assessment of the probability for a 
model is given by p = 1/(1 + B). Thus if the universe is presumed to be flat, 
the probability for the EdS model is 0.34, a reasonably good probability. In 
the absence of any prior assumption about the spatial geometry of the uni- 
verse, this probability is only 0.009 (which is still larger than 0.008 obtained 
from the frequentist approach). These results, together with the results and 
discussion of the previous section, make a reasonable case for the EdS model. 

5. EFFECTS OF WEAK LENSING 

Weak gravitational lensing is an unavoidable systematic uncertainty in the 
use of SNe la as standard candles. As the universe is inhomogeneous in 
matter distribution, the SNe fluxes are magnified by foreground galaxy ex- 
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cess and demagnified by foreground galaxy deficit, compared to a smooth 
matter distribution. Recently Williams & Song (2004) have reported such a 
correlation between the magnitudes of 55 SNe from the sample of Tonry et 
al (2003) and foreground galaxy overdensity. They have found the difference 
between the most magnified and the most demagnified SNe as about 0.3-0.4 
mag. Wang (2004) has claimed further evidence of weak lensing in the high 
redshift sample of SNe la from the Riess et al (2004) data. She has found a 
high magnification tail at the bright end of the distribution, and a demagni- 
fication shift of the peak of the distribution towards the faint end, which are 
signatures of weak lensing. She has estimated the possible magnification of 
the three intrinsically most luminous SNe (in the absence of lensing) in the 
bright end tail as: 

SN1997as (z = 0.508, p Q = 41.64): magnified by 2.10 ± 0.68; 
SN2000eg (z = 0.540, fi = 41.96): magnified by 1.80 ± 0.70; 
SN1998I (z = 0.886, /i = 42.91): magnified by 2.42 ± 1.98. 

These high magnification factors ~ 2 from weak lensing are though somewhat 
surprising, as also mentioned by Menard & Dalai (2004) who claim not to find 
any significant correlation between SN magnification and foreground galaxy 
overdensity. However, excluding these three lensed SNe from the sample 
indeed improves the fit: x7 dof = 192.49/152 = 1.27, Q = 0.015 obtained for 
M = 43.41 and np^H^ 1 = 4.82. The fit to the ACDM models also improves 
considerably: 

X 2 /dof = 164.20/151 = 1.09, Q = 0.219 obtained for fi m0 = 0.47, fi A0 = 1.02, 
with M = 43.32; 

X 2 /dof = 166.86/152 = 1.10, Q = 0.194 obtained for Q m0 = 1 - Q A0 = 0.30, 
with M = 43.34. 

Though the weak lensing effects are estimated to be small for SNe at 
z < 1, they are non-negligible for higher redshift SNe. As more SNe are 
discovered at higher redshifts, it becomes increasingly important to minimize 
the effect of weak lensing by, for example, using flux-averaging (Wang, 2000). 

6. CONCLUDING REMARKS 

Since the early 1980s, inflation and cold dark matter have been the dominant 
theoretical ideas in cosmology. However a key prediction of these ideas - the 
canonical Einstein-de Sitter model (fl m o = 1, A = 0)- seems to be in trouble 
in explaining the high redshift SNe observations, as is widely believed. How- 
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ever, it can still explain these observations successfully if one takes account 
of the extinction of light by intergalactic metallic dust ejected from the SNe 
explosions, as has been shown earlier (Vishwakarma, 2003) by considering 
Perlmutter et al' data (Perlmutter et al., 1999). We have shown, in this 
paper, that the model also has an acceptable fit to the recently published 
'gold sample' of 157 SNe la by Riess et al (2004) which, in addition to having 
previously observed SNe, also includes 16 newly discovered SNe la, 6 of them 
being among the 7 highest redshift SNe la known so far, all at redshift > 1.25. 
Though the fit to this new data is deteriorated considerably, we believe that 
one should wait for more accurate SNe la data with z significantly > 1 (as 
expected from the coming space missions such as SNAP and JWST) in order 
to rule out a model which seems to explain all other existing observations 
well, some even better than the ACDM models. This point of view is also 
corroborated by the Bayesian analysis of the data, which indicate that there 
is no significant evidence against the EdS model if one assumes a spatially 
flat universe. Even in the absence of such an assumption, there is not a very 
strong evidence against this model, for suitable priors. 

We should also examine more critically, than has been done hitherto, the 
assumption of a non-evolving standard candle for SNe la with high redshifts. 
Though most studies confirm that the luminosity properties of SNe at differ- 
ent redshift and environments are similar (Perlmutter et al., 1999; Sullivan 
et al, 2003), however, there are other theoretical studies which have found 
variations indicating evolutionary effects (Dominguez et al, 1998; Hoflich et 
al, 1998). Additionally, it has been shown by Drell et al (2000) that a com- 
parison of the peak luminosities estimated for individual SNe la by two dif- 
ferent methods are not entirely consistent with one another at high redshifts, 
z ~ 0.5. If evolution was entirely absent, the differences between them should 
not depend on redshift. Moreover the three luminosity estimators in practice 
(the multicolor light curve shape method, the template fitting method, and 
the stretch factor method) reduce the dispersion of distance moduli about 
best fit models at low redshift, but they do not at high redshift, indicating 
that the SNe may have evolved with redshift (Drell et al, 2000). 

The possible role of gravitational lensing in amplifying the supernova lu- 
minosity at high redshifts has been discussed by several authors and we have 
applied those ideas here to illustrate the difference it can make to any con- 
clusion drawn from the data. Additionally, the role of intergalactic whisker 
dust still remains to be appreciated fully and we have demonstrated here the 
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possible difference it can make to the viability of a model. 

One may think that the whisker dust can create too much optical depth 
for the high redshift objects and they need to be excessively bright in order 
to be seen. However, from our calculations, we find that even the objects 
with redshift as high as 5 will be fainter by ~4 magnitudes only, with the 
kind of whisker dust (best-fitting value) we have been talking about. Thus 
the microwave emission from the high redshift quasars will also be visible 
without unrealistically high demands on their luminosities. 

ACKNOWLEDGEMENTS 

The author thanks the Abdus Salam ICTP for sending the necessary litera- 
ture under author's associateship programme. Thanks are also due to Money 
V. John for discussion, to Victor M. Banuelos Alvarez for providing help in 
computation, and to an anonymous referee for useful comments. 

REFERENCES 

Aguire A. N., 1999, ApJ, 512, L19 

Banerjee S. K., Narlikar J. V., Wickramasinghe N. C, Hoyle F., Burbidge, 

G., 2000, ApJ, 119, 2583 
Blanchard A., Douspis M., Rowan-Robinson M., Sarkar S., 2003, A&A, 412, 

35 

Blanchard A.,2005, preprint, astro-ph/0502220 
Cayrel R., et al, 2001, Nature, 409, 691 

Chitre, S. M., Narlikar J. V., 1976, Astrophys. Space Sc., 44, 101 

Csaki C, Kaloper N., Terning J., 2002, Phys. Lett. B 535, 33 

Csaki C, Kaloper N., Terning J., 2004, preprint, astro-ph/0409596 

Dominguez I., et al, 1998, preprint, astro-ph/9809292 

Donn B., Sears G. W., 1963, Science, 140, 1208 

Drell P. S., Loredo T. J., Wasserman I., 2000, ApJ, 530, 593 

Freedman W. L. et al., 2001, ApJ., 553, 47 

Gnedin O. Y., Lahav O., Rees M. J., astro-ph/0108034 

Hoflich P., Wheeler J. C, Thielemann F. K., 1998, ApJ, 495, 617 

Hoyle F., Wickramasinghe N. C, 1988, Astrophys. Space Sc. 147, 245 

Hoyle F., Burbidge G., Narlikar J. V., 2000, A Different Approach to 

Cosmology, (Cambridge: Cambridge Univ. Press) 
Jha S., et al., 1999, ApJ. Suppl., 125, 73 
John M. V., Narlikar J. V., 2002, Phys. Rev. D, 65, 043506 



13 



Menard B., Dalai N., 2004, preprint, astro-ph/0407023 

Nabarro F. R. N., Jackson P. J., 1958, in Growth and Perfection in Crystals, 

eds. R. H. Duramus, et al, (J. Wiley, New York) 
Narlikar J. V., Wickramasinghe N. C, Sachs R., Hoyle F., 1997, Int. J. Mod. 

Phys. D, 6, 125 

Narlikar J. V., Vishwakarma, R. G., Burbidge G., 2002, PASP, 114, 1092 
Parodi B. R., et al, 2000, ApJ., 540, 634 
Perlmutter S., et al, 1999, ApJ., 517, 565 

Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 1986, 

Numerical Recipes, (Cambridge University Press) 
Riess A. G., et al., 2004, ApJ., 607, 665 
Sullivan M., et al, 2003, MNRAS, 340, 1057 
Tonry J. L., et al, 2003, preprint, astro-ph/0305008 
Vishwakarma R. G., 2002, MNRAS, 331, 776 
Vishwakarma R. G., 2003, MNRAS, 345, 545 
Wang Y., 2000, ApJ., 536, 531 
Wang Y., 2004, preprint, astro-ph/0406635 

Wickramasinghe N. C, Wallis D. H., 1996, Astrophys. Space Sc. 
240, 157 

Williams L. L. R., Song J., 2004, preprint, astro-ph/0403680 



14 



